Polarization transitions in quantum ring arrays 
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We calculate the zero temperature electrostatic properties of charged one and two dimensional 
arrays of rings, in the classical and quantum limits. Each ring is assumed to be an ideal ring of 
negligible width, with exactly one electron on the ring that interacts only with nearest neighbor 
rings. In the classical limit we find that if the electron is treated as a point particle, the ID array 
of rings can be mapped to an Ising antiferromagnet, while the 2D array groundstate is a four- fold 
degenerate "stripe" phase. In contrast, if we treat the electrical charge as a continuous fluid, the 
distribution will not spontaneously break symmetry, but will develop a charge distribution reflecting 
the symmetry of the array. In the quantum limit, the competition between the kinetic energy and 
Coulomb energy allows for a transition between unpolarized and polarized states as a function of 
the ring parameters. This allows for a new class of polarizable materials whose transitions are based 
on geometry, rather than a structural transition in a unit cell. 

PACS numbers: 73.21-b 



Advances in technology have made it possible to fab- 
ricate structures like dots, and wires on the micro- or 
nano-scale. Such quantum dots act as artificial atoms, 
with spectra and shell structure similar to those of 
real atoms. Recently it has also become possible to 
fabricate arrays of rings, each only tens of nanometers 
across. They can be fabricated either by dry etching^ 
or by using MBE techniques to foster self-assembled In- 
GaAs/GaAs rings. If a quantum dot can be thought 
of as an artificial atom, then rings (and other more com- 
plicated structures) can be thought of as "atoms" that 
do not exist in nature. We can vary their geometry to 
control properties such as polarizability, and coupling 
strength in order to produce new types of correlated sys- 
tems. 

In this letter we report the results of analytical and 
numerical calculation of polarization of one and two di- 
mensional arrays of singly charged interacting quantum 
and classical rings at zero temperature (fig.Q). The re- 
striction that each electron resides on a ring causes them 
to exhibit different ground state properties when they 
are assembled into arrays of rings than are found in ar- 
rays of simply connected quantum dots. We find that the 
minimum energy configuration of classical point charges 
in an array of rings is an antiferroelectric (AFE) pat- 
tern in ID and a striped AFE pattern in 2D square and 
triangular arrays. If we replace the point charge with 
an extended classical charge fluid, we no longer sponta- 
neously break the symmetry, but instead produce charge 
distributions that reflect the symmetry of the array. In 
the quantum limit we find that the groundstate makes 
a transition from uniform to AFE order in both ID and 
2D as a function of the array parameters. In none of our 
results were globally ferroelectric (EE) order obtained. 
We present these findings below, and then conclude with 
a discussion of how they might be tested by experiment. 

Classical Results: In order to gain intuition for the 
quantum problem, wc will examine the classical case. 



Here there are two limits: discrete charges and contin- 
uous charges We first examine the limit of ideal point 
charges on each ring with nearest neighbor interactions 
only. 

We start by considering a ID array of ideal rings. The 
radius of each ring is R, and their center-to-center sep- 
aration is D. On each ring we place one electron with 
charge e (fig.Ql). Because the electrons are confined to 
move on a ring, their moment is constrained in magni- 
tude but not in direction. 

The energy of a ID array is given by Uid = 
J2iLi ^ I l^i(^i) — where Qi is the location of 

the i-th electron as measured from the horizontal axis. 
In the dipole approximation we can write this as 
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In 



where e = RjD and C/g is a constant, \Jq - 

the second expression we write the position of the charge 
as Si, a vector in the 2D plane pointing from the center 
of the z-th ring to the charge on that ring. The cos 2Q 
(or D ■ s) term explicitly breaks the rotational symme- 
try, driving the system from XY to Ising-like behavior. 
The last two terms in eq.(^ compete: the first driving 
the system to be ferroelectric (like direction) while the 
latter, larger term favoring states where neighbors point 
in opposite directions. Thus the system at zero temper- 
ature will order in an AFE pattern in ID. This is born 
out by numerical Monte Carlo simulations of the exact 
Coulomb interaction, which reproduce the AFE pattern 
on finite size arrays (Fig^J. 

If we expand the potential to quadratic order in the 
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spectrum is given by 



FIG. 1: Figure 1: A scliematic picture of the groundstate of 
classical point electrons for (a) ID and (b) 2D square arrays of 
rings. The ring radius is 7? and the separation is D. The ID 
ordering is antiferroelectric and thus has a double degenerate 
groundstate. The 2D square array is striped antiferroelectric, 
and thus has a fourfold degenerate groundstate. 



displacement and calculate the normal modes we find: 
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where c^o = ^Je^ /m*D'^. Both the modes are gapped, 
since there is an Ising-likc term that provides a harmonic 
restoring force at each site. Note that they are indepen- 
dent of ring radius at this level of approximation.. 

We may repeat the analysis above for a 2D array of 
rings. The Coulomb energy for an N x N array may be 
written as 
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where i and j denote the x and y indices of the ring in 
the plane. Since each ring in the plane (except the ones 
on the top and right boundaries) has one neighbor to 
the right and one above it, this covers the array without 
double counting. If we neglect edge effects the minimum 
energy configuration of the system is AFE in one direc- 
tion (x or y) and FE in the other (y or x). Since the 
AFE state is itself doubly degenerate, we have a four- 
fold degenerate minimum energy state (Fig.(Q. This is 
supported by numerical simulations which enforce (un- 
physical) periodic boundary conditions on the array. In 
more realistic 2D finite arrays the long range order is 
perturbed by a "vortex" like pattern on the boundaries 
where the charge on each ring resides on the point that 
puts it on the outermost perimeter of the array. 

A normal mode analysis shows that the excitation 
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where we have assumed that the stripes run vertically in 
the ordered phase, as in fig.lQJ. Note that the acoustic 
mode is now ungapped, at the wavevector fco = {0,tt/D), 
in contrast to the ID case. This result holds in the case 
of only nearest neighbor interactions, where the increase 
in the energy due to departure from AFE order in the 
horizontal neighbors is compensated by the motion of 
vertical neighbors from FE allignment. This mode differs 
from the magnetoplasmonQ mode of metal ring arrays. 

A second classical limit that contrasts with the point 
charge model is a continuous charge fluid. That is, 
each point on the ring has a charge density pi{9i) such 
that the integrated charge density per ring is a constant, 
/ Pi{6i)d6i = 1. We seek a minimal solution to the vari- 
ational quantity: 
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For a ID ring this is expression is divergent due to the in- 
finite self energy. We can regularize this in several ways. 
One method is to introduce a short distance cutoff A to 
the Coulomb interaction, discretize the integral equation, 
and then solve the problem numerically. This can also be 
done analytically, expanding all results to lowest order in 
the cutoff, e. 

An approximate analytic solution can then be obtained 
by Fourier expanding the distribution, keeping only the 
first three modes. We find that amplitude non-trivial 
Fourier mode as a function of A and e = R/D for p is 
given by 
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We compare this analytic result with the numerical ap- 
proach in Fig.||2Jl. A second method is to obtain analytic 
results by smearing the ring charge density over a very 
short cylinder, whose height h is small compared to its 
radius. Analytic results are in qualitative agreement with 
the cutoff expansion. 

Quantum mechanical results: The quantum mechani- 
cal case at first glance seems similar to the classical fluid 
charge case. Although there is no self-interaction, the 
quantum particle has a kinetic energy which opposes lo- 
calization, hence for low values of interaction strength 
quantum fluctuations should destroy the polarization. 
However for large values of interaction it is not a pri- 
ori obvious whether the polarization pattern beaks the 
symmetry. 

Again we assume exactly one charge on each ring and 
only nearest neighbor Coulomb interaction. The thick- 
ness of each ring is much smaller than its radius so that 
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FIG. 2: A plot of the second fourier amplitude of the classical 
charge distribution on a ring in a ID horizontal array. The cir- 
cles are numerical results, the solid line is a scaled plot of eq|K| 
Scaling is required since the analytic result neglects all higher 
fourier modes. Inset: a sketch of the charge distribution that 
corresponds to this Fourier mode. Note that the symmetry of 
the array is not broken by the charge distribution. 

the only coordinate of the charge is the angular position. 
We decompose the wavefunction in each ring into a lim- 
ited number of Fourier modes, tpi{0) — J2n=~no ^^e^^^ , 
and then solve the system numerically in the Hartree ap- 
proximation. We use a dimensionless Hamiltonian for the 
i-th ring: 



de' (6) 



where the sum is over the nearest neighbors of i-th ring 
and 5 = (e^ /D) / (fi^ /2mR^) is the interaction strength 
and energy is measured in units of /2mR^. We im- 
pose periodic boundary conditions on the array and by 
an iterative self-consistent method we find that in a ID 
array of rings for a finite value of 6 — Sc there is a change 
of ground state from totally unpolarized state |P| = 
to AFE state |P| = Pq where P is the staggered polar- 
ization vector P = ^Ei(-l)'/ ni9)\i:,{e)\'^d9 for ID 
array. The ground state energy per ring: 
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also shows a change in behavior at 5c- Results of these 
calculations are plotted in Fig.©. 

We can obtain analytical insight by assuming that the 
ground state wave function for all rings in each sub- 
lattice is the same, so that we can employ a simple 
variational wave function for each sublattice ^a(^) = 



^cos(6'-0) and V^b(6') 



alternatively. Inserting these functions into (0 and min- 
imizing the energy using dipole expansion we find an ap- 
proximate behavior for energy and polarization. In the 




FIG. 3: A comparison numerical and analytical calcula- 
tions of the staggered polarization as a function of 5 = 
in a ID quantum ring array obtained 
in the Hartree approximation. The numerical results are for 
the case include Fourier modes |m| < 1 (dashed line) and 
|m| < 6 (circles). The solid line is the analytic result assum- 
ing |m| < 1. The quantity 5 is a measure of the competition 
between the Coulomb interaction and the quantum kinetic 
energy. 



ground state we find 4> = ^ and 




for 5 > Scie) = ^e-^ 
for 5 <5c 
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The numerical Hartree results for the energy and po- 
larization are in agreement with the variational calcula- 
tion when we restrict the number of Fourier modes in 
the numerical calculation to n € { — 1,0,1}. The results 
are little changed when more Fourier modes are intro- 
duced to the numerical calculation or even when the ex- 
act Coulomb interaction is used in our numerical code. 

Using the same Hamiltonian as eqn© for a 2D array 
and with periodic boundary conditions we are able to find 
the behavior of the ground state energy and staggered 
polarization vector as a function of interaction strength. 
Again by using limited number of Fourier modes for each 
wave function we find that there is a change of ground 
state from totally unpolarized state to a striped AFE 
state in which each column has FE polarization but in 
opposite direction of its neighbor column. The staggered 
polarization vector for this veritically oriented case is 
expressed as .P = J2i ji~^yPij/^ where j counts the 
columns and pij = J rij(0)\ {9)\'^d9 is polarization vec- 
tor of each ring. The ground state is four-fold degenerate 
(the states obtained by flipping every spin and or rotat- 
ing them all by 7r/2 yield the same energy) and the U(l) 
symmetry of unpolarized state has been broken. 

Again it is possible to estimate a variational wave func- 
tion for each ring assuming the wave functions of the rings 
in each sublattice are the same (infinite system or peri- 
odic boundary condition). Using dipole approximation 
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FIG. 4: A comparison numerical and analytical calcula- 
tions of the staggered polarization as a function of 5 — 
(e^/D) / (h^ /2mR^^ in a 2D quantum ring array obtained 
in the Hartree approximation. The numerical results are for 
the case include Fourier modes |m| < 1 (squares) and |m| < 6; 
the dashed lines are added as a guide. The solid line is the 
analytic result assuming \m\ < 1. The quantity (5 is a measure 
of the competition between the Coulomb interaction and the 
quantum kinetic energy. 

we find = f and 

[ 0, for 6<Sc 

Numerical Hartree data and variational approximation 
results are in good agreement even by using exact 
Coulomb interaction in numerical calculations FigQ). 

The quantum phase transition in the ID system should 
be similar to that of the transverse field Ising model. Q 
The cos 26 term breaks the XY-symmetry so that at low 
temperatures the system is Ising. The kinetic energy 
term tries to delocalize the spin, similar to the transverse 
magnetic field in an Ising model. However, the Coulomb 
interaction is not precisely an AFE Ising interaction; it 
is a difference of cos {6 — 6') and cos {6 + 0') terms that 
we can cast it in that form only at very low temperatures 
where we assume that < 9 > \s sharp on each ring. 

In 2D there should be a finite temperature phase 
transition between the polarized and unpolarized state. 
While this Letter focuses only on the zero temperature 
behavior, Monte Carlo simulations of the I+ID quantum 
problem not only support the existence of a ID quantum 
phase transition, they document the existence of a 2D 
transition at finite temperature. 

Considerations for experiment: The parameter 5 de- 
termines if the array will spontaneously polarize; in ID 
the transition is at (5 w 10, while in 2D, it is at (5 w 5. 
It is easy to achieve small values of 5 simply by choos- 
ing the ring separation to be large. Thus the "quan- 
tum" limit where the kinetic energy dominates is sim- 
ple to obtain. To obtain the antiferroelectrically or- 
dered state we need large 5. We may write this as 



5 = [R^/aoD) X {m* /m) where oq is the Bohr radius 
and m* is the effective mass of the electron. If R and D 
are R and D measured in nanometers, and m = {m* /m), 
then 6 « 18. 9 {R^ / D)rh. We require that the rings do not 
intersect, so that D > 2R. Thus the ability to achieve 
large values of S in semiconductors will depend upon the 
value of the effective mass. If we set D = 2R, then for 
GaAs (m — 0.06) ID arrays of rings with a radius greater 
than ~10nm will be polarized. For AlAs (m — 0.4) the 
crossover radius is about 70nm. Rings with a smaller 
radius will not spontaneously polarize, but instead be 
isotropic 

It is well known that in ID there is no ordered state 
for r > for the Ising model. However, this is true 
in the thermodynamic limit. For small arrays over fi- 
nite time intervals the system can order. If we wish 
to observe this behavior we want the characteristic en- 
ergies of the system to be greater than the tempera- 
ture. For the Coulomb energy kT < e^/D, which we 
may write as DT < 1.8 x 10"^ were D is measured in 
nm and T in Kelvin. For the kinetic energy this means 
kT < h/2m*R^; if we measure {m/m*)R^T > 40 in the 
same units. For GaAs we can choose R to be about 14nm 
at 4K; choosing materials with a smaller effective mass 
or going to lower temperature allows us to increase the 
radius. 

A AFE polarized ring array will scatter light at a 
wavelength commensurate with the inter-ring separation, 
D. In ID there is a gap 1/2^0, which we may write 
as 2^/2rh{ao/D)^^^. For GaAs rings with a separation 
D = lOOOnm this gives - 6.0 x 10^°Hz. The 2D arrays 
have a similar sized gap at zone center, but the gap van- 
ishes at one zone edge. The excitation spectrum can be 
probed optically, but scattering at the edge of the zone is 
difficult due to the constraints imposed by conservation 
of energy and momentum. Typically in such cases Ra- 
man scattering can be used to investigate the excitations. 

Finaly these calculations assume that each ring is 
singly occupied. This might be obtained by fabricating 
the rings upon a thin insulating layer covering a gate. By 
tuning the gate voltage we can bias the system so that it 
is energetically favorable for an electron to tunnel to the 
rings. The gate will also serve to cutoff long distance in- 
teractions between the rings, supporting the assumption 
of the nearest neighbor interactions used here. Moreover, 
this letter serves to start investigation into a broad class 
of problems, such as rings occupied by an optically ex- 
cited exciton/hole pair or perhaps by a small, varying 
number of electrons created by a random distribution of 
dopants. 

The topic of quantum dot arrays and their correlations 
has obvious and useful analogies with solid state models 
of crystalline arrays of atoms. In this Letter we wish to 
point out that experimentalists have at their disposal a 
host of "unnatural atoms" analogs: rings, quantum dot 
quantum wells, quantum rice, etc. The electrons in these 
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nanoscale constituents are confined to orbitals that may 
not have atomic analogs. Morever, it may be possible to 
tune the shape of the constiutent to optimize some de- 
sired collective property such as frustration in electric or 
magnetic polarization, high susceptibility or sensitivity 
to optical polarization of light. Even more rich behavior 
will develop if we allow electrons to tunnel between these 
nanoscale periodic structures. 

In this Letter we have examined the simple case of 
arrays of singly charged classical and quantum rings at 
zero temperature. We have found AFE order in ID clas- 
sical systems and a transition to ID AFE order in the 
ID quantum systems. In 2D classical and quantum sys- 
tems there is a similar striped ordering. The ordered 
phases have collective excitations that can be measured 
optically. 
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